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Abstract 

We construct a mixed spin 1/2 and S integrable model and investigate its finite 
size properties. For a certain conformal invariant mixed spin system the central 
charge can be decomposed in terms of the conformal anomaly of two single integrable 
models of spin 1/2 and spin [S — 1/2). We also compute the ground state energy 
and the sound velocity in the thermodynamic limit. 
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Integrable magnetic spin chains provide important examples of systems which can be 
derived from the so-called Yang-Baxter algebra |Q. A well known model is the isotropic 
spin 1/2 Heisenberg ^ chain and its generalization for arbitrary spin 5 Another 
interesting example is the Heisenberg model in presence of an impurity of spin 5 |^. 
In such model one of the local vertex weight acts on a pair of asymmetric vector spaces 
which is defined by the local states on the horizontal and vertical lines of a two dimensional 
lattice. More recently, a general discussion concerning the construction of mixed vertex 
models has been presented by de Vega and Woynorovich [0]. For instance, they have 
studied several properties of the thermodynamic limit of an alternating anisotropic chain 
of spins 1/2 and 1. However, it is still to be investigated the finite size effects in these 
mixed spins models as well as their class of universality for conformally invariant systems. 
Following the approach of ref. |^ we construct an isotropic alternating spin 1/2 and S 
chain. We focus our attention in the analysis of the finite size behaviour of the ground 
state on a line of length L. A conformally invariant mixed system can be defined and 
its conformal anomaly is computed by analysing the finite size corrections for the ground 
state energy and by the thermodynamic Bethe ansatz. Several useful quantities as the 
ground state energy and the sound velocity are also computed. 

The construction of the transfer matrix of the mixed spins a and 5* model is based 
on the local vertex Rgj{X) which is a matrix in the auxiliary space Va and its matrix 
elements are operators of spin 5* acting on the Hilbert state space at site j. In the case of 

1 /2 

an auxiliary space of spin 1/2, Rgj (A) is given by 

1/2 ( A^o + iS^ iS~ \ _^ , 

RsiW = Sf = S]±zS^ (1) 

where Sj, a = 1, 2, 3 are spin 5* operators and 5*° is a (25 + 1) x (25* + 1) identity matrix. 

In the case of an alternating spin 1/2 and 5* mixed model the set of commuting transfer 
matrix Ti/2,s(A) assumes the following form 

Ti/2,5(A) = Trv.,,,(n/2,5(A)), n/,,sW = JA)<Li(A) ■ • ■<^(A) (2) 



E E <-i{si si}cri^, + E(^ - s{s + 

n=even i,j=0 i=l 



where Ti/2,5 is a (2x2) matrix in the auxihary space V1/2 denominated the monodromy 
matrix. Here we impose periodic boundary conditions and the length L is an even number. 

The associated one dimensional quantum Hamiltonian which commutes with Ti/2,5 is 
defined by i/1/2,5 = log(Ti/2,5(A))|A=i/2 and has the following expression 

JL{2S + 3)^ 
8 

(3) 

where a*, i = 1, 2, 3 are Pauli matrix elements, cr° is the identity matrix, and J = 23/ (25* + 
1)^ . In this paper we are interested in the antiferromagnetic (J > 0) regime of (^) and 
let us assume, for the sake of simplicity, J = 1. 

Similarly to the usual spin 1/2 Heisenberg chain, the Hamiltonian can be diagonal- 
ized by the quantum inverse scattering method 0. The eigenenergies are parametrized in 
terms of the complex numbers \j 



M 



where \j satisfy the so-called Bethe ansatz equation 

where M = ^{S + 1/2) —r, and r labels the disjoint sectors of the eigenvalues of the total 
spin operator Ejlodd(^i) + T.f=even{S j) ■ 

In order to find the structure of the numbers \j for finite size systems we numerically 
solve Eqs.(4,5) and compare them with the exact diagonalization of the Hamiltonian (3). 
In Fig.l(a,b) we show the picture of the ground state for L = 6, 8, 10 and S = 1, 3/2. In 
Table 1 we present the respective values of the complex parameter Aj. For large enough 
L (L > 8) the numbers Aj cluster in two distint sets of roots. One of them consists of 
real numbers (full circles) and the other of complex structures (crosses) A" which in the 
asymptotic limit L 00 are called 2S-strings 

AJ = e, + ^(2S + l-2«), « = 1,2,...,2^ (6) 



where is a real number denominated center of the 2S-string. 

Taking into account the structure discussed above, the ground state energy per particle 
e^^'"^ can be calculated and is given by 

eV2- = - ln(2) - 1 (^(^) - ^(^)) (7) 

where iIj{x) is the Euler psi function. 

Analogously, it is possible to define the transfer matrix T5' i/2(A) of a mixed spin S and 
1/2 model which commutes with Ti/2,5(A). The transfer matrix Ts^i/2{^) has the following 
expression 

T5,i/2(A) = Trvsirs,,/2iX)), rs,i/2{X) = i?f/2, JA)i?f,^_i(A)....i?|,(A) (8) 

The commutativity between Ti/2,5(A) and T5' i/2(A) derives from the Yang-Baxter rela- 
tions satisfied by the local vertices Rgj 0. Interesting enough, this permits us to define a 
rotational invariant mixed vertex model by formally multiplying the two transfer matrix 

i 

T^ym = T^/2,sTs,i/2 (9) 

Due to commutativity [Ti/2,5(A), 1/2A'] = 0, the eigenenergies of the one-dimensional 
Hamiltonian associated to T'^^™' are parametrized by the same Bethe equations, namely 
Eq.^ However, the expression for the total energy for a given set of numbers Xj is 

M -I M cyn 

^•'--E^-E^ (10) 

The model defined by Eqs.(|, [1^) possesses all features of being conformally invariant, 
i.e, short-range interaction, translational and rotational symmetry and gapless low-lying 
excitations in the spectrum. Indeed for low (total) momenta p Q,the dispersion relation is 
linear in p 

e{p) ~ VsP (11) 
^We remark here that the correct momentum operator is half of that used in ref.[7] (see eq.(2.35)). 
This imphcs that the sound velocity is double {vg — 2tt) of that found previously in ref. [7] (vs — tt) . 
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where Vg = 27r is the sound velocity. 

The class of universality of this model can be determined by exploiting a set of impor- 
tant relations between the eigenspectrum of finite lattice systems. In particular the 
conformal anomaly is related to the ground state energy Eq^"^{L) by ||Tl], |12[ 



EriL)/L = - ^ (12) 



where c is the central charge and e^™ is the ground state per particle 

1 



sym _ 1/2,5 



,,25 + 3, ,,25+L ,,25 + 1, , / 

^i—^) - H—^) + H^^) - V'(l/2) 



(13) 



In Table 2 we present our estimatives for the central charge c of Eq.(|T2|) for S = 
1,3/2. Our numerical result predicts a conformal anomaly c = 2.01(2)(S' = 1) and c = 
2.500(6)(S' = 3/2). We notice that for 5 = 1 the extrapolation is less precise compared to 
S = 3/2. This is due to the fact that the bulk of the complex part of 2-string zeros are 
next to ±i/2. As a consequence, we can apply an analytical technique developed by de 
Vega and Woynorovich |]13| and we find the exact value c = 2 for S = 1. 

Another efficient method to compute the conformal anomaly is by analysing the low 
temperature behaviour of the associated free energy. For a critical system at low temper- 



ature, the free energy per particle has the following asymptotic behaviour |TT], [12 



F{T)/L = eoo - -— (14) 

In the case of integrable one dimensional spin chains the thermodynamic properties 
can be studied by using the thermodynamic Bethe ansatz method. In this approach, the 
free energy is given in terms of variables denominated pseudoenergies and its minimization 
yields a set of integral equations for these parameters. In order to obtain such equations 
for the mixed spin model defined by Eqs.(|^, |iy) we follow refs. p, Q and here we give 
only the final results. The free energy at temperature T is given by 

F'y"'{T)/L = e^™ - — / p{X) (ln(l + e^^W/^) + ln(l + e'^^^^^^^) (15) 

4:71 J-OO ^ ' 
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where the pseudoenergies ea(A), a = 1, 2, 3, . . . satisfy the following TBA equations 

en{\)/T = p{X) * (ln(l + e-+^/^) + ln(l + e--/^)) + ^(5„,i + 6n,2s) (16) 
where eo(A) = 0, f*g{x) denotes the convolution ^ f{x)g{x—y)dy andp(A) 



cosh(7rA) ■ 



The advantage of Eqs.(|l5|, p^ ) is that they allow us to exactly estimate the low tem- 



perature behaviour of the free energy. Using the standard procedure ( see e.g. [15, ffl) to 



compute the leading behaviour of F{T), and some dilogarithm identities we find the 
following result 



Comparing Eqs.(|T7D and ([T^) the value of the central charge is given by 

For S = 1, 3/2 we obtain c = 2, 2.5 which are consistent with our numerical findings of 
Table 2. Interesting enough, Eq.(|18D can be decomposed in terms of the central charges 



of the integrable Heisenberg chains of spins 1/2 (c = 1) and S - 1/2 (c = ^g^). In 
this sense the effect of the interaction between spins 1/2 and S is the "reduction" of the 
critical behaviour of a spin S" to S* — 1/2 0. Taking into account this last observation it is 
easy to conjecture the conformal anomaly of a mixed spin S and S' . In general, we can 
assume S' > S and the expected critical behaviour is given by the following conformal 
anomaly 

3(S'-S) 

S + l^ S'-S + l ^ ' 

We believe that result (|19|) is the first step toward the understanding of the composition 
of the operator content of these mixed spins models. However, it is still to be investigated 
how the primary fields are composed in these systems. We hope to report on this problem 
in a future publication. 



^ We remark that J = S"— 1/2 is the smallest possible value in the addition of two "angular momenta" 
S and 1/2. 
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Figure Captions 

Fig.l(a,b) The ground state configurations of the numbers Xj for L — 6,8,10 and 
5 = 1,3/2. 



Table Captions 



Table 1. The numerical values of the numbers Xj of the Bethe ansatz equation (5) corre- 
sponding to the ground state for L = 6, 8, 10 and S — 1, 3/2. 

Table 2. The estimatives of the conformal anomaly of Eq.(12) for 5* = 1 and S — 3/2 
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Table 1 



L 


S = l 


5 = 3/2 


6 


0.316194 ±i0.508074 
-0.494077, 0.228965 


±il.5, ±i0.5 
±0.358639 


8 


±0.308362 ± i0.503329 

±0.308005 


±0.302413 ±il. 060323 

±0.436072, ±0.201077 


10 


±0.535111 ±i0.510182 
±0.444828, 


±0.302415 ±il.060320 
±0.436069, ±i0.5 
±0.201079 



Table 2 



L 


5 = 1 


5 = 3/2 


8 


2.210 923 


2.839 364 


16 


2.061 249 


2.602 543 


24 


2.041 930 


2.556 301 


32 


2.020 828 


2.538 530 


40 


2.017 593 


2.529 503 


48 


2.014 692 


2.524 142 


56 


2.012 134 


2.520 624 


Extrapolated 


2.01(1) 


2.500(6) 
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(a) 



L = 6 



L = 8 



L = 10 




Figure l(a,b) 
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